Tunneling density of states of granular metals 
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We investigate the effect of Coulomb interactions on tlie tunneling density of states (DOS) of 
granular metallic systems at the onset of Coulomb blockade regime in two and three dimensions. 
Using the renormalization group technique we derive the analytical expressions for the DOS as a 
function of temperature T and energy e. We show that samples with the bare intergranular tunneling 
conductance g% less than the critical value g^ — {1/2-Kd) ln{Ec/5), where Ec and 5 are the charging 
energy and the mean energy level spacing in a single grain respectively, are insulators with a hard gap 
in the DOS at temperatures T — > 0. In 3d systems the critical conductance g^ separates insulating 
and metallic phases at zero temperature, whereas in the granular films g^ separates insulating states 
with the hard (at g^ < g^) and soft (at g% > gr ) gaps. The gap in the DOS begins to develop at 
temperatures T* ~ EcgT exp{—2TTdg^) and reaches the value A ~ T* at T — » 0. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.+h 



Granular materials, a focus of the current mesoscopic 
physics, is a unique testing, area for general concepts 
of disordered systems (see S of the re- 

markable fundamental features of granular metals is the 
strong on-site Coulomb repulsion that leads to the sup- 
pression of transport at low temperatures and to an in- 
sulating ground state at low bare intergranular conduc- 
tances < 5t = (l/27rd) ln(i?c/<5). A good progress 
in understanding this insulating state has recently been 
made Q; yet, despite the recent advance, the satisfactory 
picture of physics of and near the Mott transition that 
occurs at — 9ti "^^^ still lacking. One of the funda- 
mental questions that remained open was the suppression 
of the tunneling density of sates that always accompanies 
Mott transition. In this Letter we develop a quantitative 
approach that enables us to investigate Mott transition 
in granular metals and derive the associated tunneling 
density of states in its vicinity. 

The density of states is a fundamental quantity that 
determines most of the properties of the system involved, 
and the electronic transport is a key phenomenon where 
the manifestations of its features may be most pro- 
nounced. A general technique to treat transport prop- 
erties of granular metals in the high temperature regime 
T > grpS was developed recently in Ref. 0. It was, in 
particular, shown that the conductivity of granular met- 
als can be described in terms of the renormalized tem- 
perature dependent intergranular tunneling conductance 
given by the following expression 



g,(T)=5° -(l/27rd) HglEa/T], 



(1) 



which holds as long as griT) > 1. The conductivity of 
the sample is related to the tunneling conductance as 
a{T) = 2e^gT{T)a?~'^, where a is the granule size, d is 
the dimensionality of the granular array and the factor 
of 2 is due to the spin. From Eq. |^ follows that at 



temperatures 



Ecgl e 



(2) 



the renormalized conductance, griT), is strongly sup- 
pressed and approaches small values where renormaliza- 
tion group breaks down. Equation is valid only at 
temperatures T > g^S whereas in the opposite case, 
T < (7^,(5, the conductance renormalization ^ is sat- 
urated and the system behaves essentially as a homo- 
geneous disordered metal 0. Comparing two relevant 
energy scales T* and grS, one concludes 3 that if (i) 
T* < grS (or, equivalently, g° > g'^) then, the renormal- 
ized conductance is still large at temperatures T ^ g^ 5 
and the low temperature phase of the system is similar 
to that of the disordered metals. Alternatively, if (ii) 
T* > g^S (or ,g° < g^), the conductance of the system 
becomes significantly suppressed at T T* reflecting 
thus the onset of the Coulomb blockade regime. In the 
latter case one expects that at T ^ T* the Coulomb gap 
begins to develop (reaching its maximal value at zero 
temperature) and, as a result, a noticeable suppression 
of DOS even at finite temperatures, T ^ T* occurs. 

In this Letter we consider the tunneling DOS of granu- 
lar metals with the bare tunneling conductance (?° < g'^ 
at the onset of Coulomb blockade regime at temperatures 
T > T* [case (ii) above] . We show that Coulomb block- 
ade strongly suppresses the tunnehng DOS, J^(T). For 3d 
granular samples i^(T) is given by 
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whereas for granular films we obtain 
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(3a) 



(3b) 



Here Vo is the DOS for non interacting electrons and 
A = 0.253 is the dimcnsionlcss constant. Equations Q 



2 



hold for temperatures T > ma.x(T* , g^S) where the tem- 
perature T* is given by Eq. (0) ■ We note that according 
to Eqs. |(2J) the DOS vanishes exactly at the same tem- 
perature T* as the renormalized conductance gxiT) in 
Eq. CEJ. 

We show that the results © can be generalized to 
finite frequency by substitution T —^ max{T, e}. In this 
case Eqs. © can be applied even for T ^ provided e 
is larger than the characteristic energy A that coincides 
with the temperature T* 



(4) 



From Eqs. JSJ one can see that the tunneling DOS 
is strongly suppressed at energies e ~ A. Thus, one 
concludes that for < g^ the system is an insulator 
at zero temperature with the Coulomb gap A given by 
Eq. (0J), see Fig. 1 a). On the contrary, for e = but at 





FIG. 1: Schematic behavior of the Density of States as a 
function of a) energy and b) temperature. 

finite temperatures, T we expect the exponentially small 
but finite activation contribution for DOS 



-A/T 



(5) 



see Fig. 1 b). 

Equation (1) generalizes straightforwardly to finite fre- 
quencies by the substitution T — > max{T, w}; this allows 
us to relate the frequency dependent conductivity, cr(w) 
with the tunneling density of sates at zero temperatures, 
T = 0: For 3d granular samples we obtain the following 
scaling relation 



3A 



(6a) 



where ag — le^g^a} is the high temperature conduc- 
tivity. For granular films we get the following expression 



ea{uj) 



(6b) 



where e « 2.7182. Equations © are useful for the com- 
parison of our predictions with the experimental data. 



Now we turn to the quantitative description of our 
model and derivation of Eq. Q: We consider a 
d— dimensional array of metallic grains. The motion of 
electrons inside the grains is diffusive and they can tun- 
nel between grains. We assume that in the absence of the 
Coulomb interaction, the sample would be a good metal. 

The system of weakly coupled metallic grains is de- 
scribed by the Hamiltonian 



H — Hq + He + Ht 



(7a) 



The term Hq in Eq. H7a(l is the Hamiltonian of an isolated 
disordered granule. The term He describes the Coulomb 
repulsion and is given by 



Hr = 



(7b) 



where i stands for the granule number, is the capac- 
itance matrix and hi is the operator of electrons number 
in the i-th granule. The last term on the right hand side 
of Eq. H7a|l is the tunneling Hamiltonian 



(7c) 



where a| fe(ai,fc) are the creation (annihilation) operators 
for an electron in the state k of the i-th grain and t.y is 
the tunneling matrix element corresponding to the points 
of contact of i-th and j—th granules. 

As it was shown in Ref. |5|, at temperatures T > gj-S 
the model given by Eq. (fT)) can be effectively described in 
terms of the functional proposed by Ambegaokar, Eckern 
and Schon (AES) in Ref. 



S ^ Sc + St 



where Sc is the charging energy 



dT- 



i4)i Cij d(j)j 
d7 2e2 dr ' 



(8a) 



(8b) 



and the second term in the right hand side of Eq. isal 
describes tunneling between the granules 



St = 2TrgT 



T^drdr' 



sin^[7rr(T-T')] 



■ sm 



'/'»j(t) - </',j(t') 



(8c) 

Here (pijir) = <pi{T) ~ (j)j{T) is the phase difference be- 
tween the phase in the i-th and j-th granules. In the 
metallic regime we may neglect winding numbers in the 
phases. We use the renormalization group technique to 
calculate DOS. The charging part, Sc of the action in 
Eq. H8a|l determines the upper frequency cutoff. 

In terms of the functional approach based on the action 
l|5a|l . DOS is expressed as (J| 



1^0 



= Tim 



dre'' 



sin ttTt 



-i(0.(r)-0,(O) 



(9) 
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Following the standard renormalization group procedure 
we separate the field 4> into slow, and fast, parts 
and integrate over the fast field in Eq. @ . The fast fields 
belong to the interval dS that represents the energy shell 
dA in the configuration space of the quasimomentum, q 
and frequency, uj. Using the one- loop approximation we 
obtain the following RG equation for the conductance gr 



dgriA) 



1 



(10) 



dA 2ndA ' 

Performing the integration in Eq. Hl()|l we obtain Eq. 1^. 
With the same accuracy for the flow equation of the den- 
sity of sates we get 



d ln(z^/z^o) = o," 



duo d'^q 



dS 2vr (27r) 



^G4^,q). (11) 



Here the Green function G^{uj, q) of the phase fields (p is 
defined on the scales A 

1 1 



Gct,{uJ,q) 



(12) 



2.gT(A) \lu\E,- 

In equation l|12|) we introduced the notation Eq = 
2 X^al-i ~ cos(qa)]. The integration in Eq. ffTT^ is going 
over the infinitesimal volume dS in the (w,q) configu- 
ration space that corresponds to the energy interval dA. 
The proper way to chose a particular form of dS depends 
on the dimensionality of the sample: For Sd samples the 
integrals over the quasimomentum converge and one can 
simply choose dS = (2Tr / a)^ dA. This leads to the follow- 
ing differential equation 

, / \ ^ dA / 1 o \ 

dln(./.o) = ^^^^, (13) 

where A — a'^ f -M-tt-^ ~ 0.253 is the numerical con- 
stant. Integrating over the A in Eq. I|13() in the range 
(T, g° Ec) we obtain Eq. l|Sa|l for the DOS of 3d granular 
metals. 

The 2d case is different since the direct integration over 
the quasi-momentum, q in Eq. (|ll|l would lead to the 
infra-red divergence. In this case it is natural to intro- 
duce the infinitesimal volume, dS in the following way 



dujd q 



dS 



dujd'^qS{\uj\Eg ~ A)dA, (14) 



such that on the energy shell A the propagator H12|l will 
not be divergent since ujEq — A. Using Eq. H14|) . per- 
forming the integration over the uj, q in Eq. (|ll|l and 
taking into account the fact that the upper cutoff for 
frequency, to is grEciq), where for 2d granular samples 
the charging energy, Ec{q) is given by the expression 
Ec{q) = nEc/qa, with the logarithmic accuracy we ob- 
tain the following equation 



dln(zy/i'o) = 



1 



1 



47r2 A5r(A) 
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9tEc 



A 



dA. 



(15) 



Integrating Eq. H15() over the variable A in the interval 
{T,gT.Ec) we obtain DOS as given by Eq. (|3b|) for 2d 
granular samples. Equation Ij3b|l was derived for long 
range Coulomb interaction. For short range Coulomb 
interaction (this may happen due to the presence of ex- 
ternal screening) the upper cutoff for w-integration in 
Eq. 1)111) is q— independent and is given by grEc. In 
this case the extra factor 1/2 will appear in the right 
hand side of Eq. 1)15)1 and as a consequence the fi- 
nal result for the DOS in Eq. ()3b)l will be modified as 
v{uj)/i'a [i>{u})/voY/'^. In this case for > g^ one 
reproduces the result of Ref. Q for the tunneling DOS 
of granular metals in the limit of large tunneling conduc- 
tance. 

Although our theory applies to Id granular arrays, the 
results in this case should be taken with caution when 
compared with experimental data: the problem is that 
the conductivity of Id system is usually controlled by 
the weakest junction and thus must be described by the 
conductance distribution function rather than by the 
average value of the conductance. Thus even small fluc- 
tuations of conductance could be important in Id case, 
especially close to the predicted transition at 5° ~ g^; 
we, however leave the detailed analysis of this situation 
to the forthcoming publication. 

In conclusion, we have investigated the effect of 
Coulomb blockade on the tunneling DOS of granular 
metals in the limit of large tunneling conductance be- 
tween the grains. We have determined the critical value 
of tunneling conductance = {l/2'Kd)\w[Ec / 5) below 
which the granular metal becomes an insulator with a 
"hard" gap at zero temperature. For 3d samples this 
value of critical conductance corresponds to a metal- 
insulator transition, as granular sarnples with g° > g'^ 
are metallic at zero temperature This value of 

= (l/27rd) ln(ii^c'/^) explains the long known puzzling 
fact that in 3d systems the metal insulator transition oc- 
curs at 5^ « 0.1. 

The situation is different for 2d granular systems since 
in this case even samples with g'^ > g^ arc insulators 
at temperatures T — > due to interaction and quantum 
effects 6], similar to those that take place in homoge- 
neously disordered metals J, 10, JTJ. Nevertheless even 
in 2d case the critical value of conductance g^ represents 
the boundary between two physically different regimes at 
temperatures T — > 0: Samples with < represent 
the "hard" insulators, with a hard gap in the DOS, while 
samples with > arc insulators with a soft gap in 
the DOS similar to homogeneously disordered metals. 
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